function b = soll(L,b)
n = length(b);
b(1) = b(1)/L(1,1);
for i=2:n
  b(i) = (b(i) - L(i,1:i-1)*b(1:i-1))/L(i,i);
end

